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The safety analysis code SCTRAN for SCWR (Super Critical Water Reactor) is modified to own the capability 
to assess the radiation heat transfer with developing a two-dimensional heat conduction solution scheme and 
incorporating a radiation heat transfer model. The verification of the developed radiation heat transfer model 
is conducted through code-to-code comparison with CATHENA. The results show that the modified SCTRAN 
code is successful for that the maximum absolute error and relative error of the surface temperature between 
results of SCTRAN and CATHENA are 6.1 °C and 0.9%, which are acceptable in temperature prediction. Then, 
with the modified SCTRAN code, the loss of coolant accident with a total loss of emergency core cooling system 
(LOCA/LOECC) of Canadian-SCWR is carried out to evaluate its “no-core-melt” concept. The following 
conclusions are achieved: 1) in the process of LOCA, the decay heat can be totally removed by the radiation 
heat transfer and the natural convection of the high-temperature coolant, even without an intervention of ECCS 
(Emergency Core Cooling System); 2) The peak cladding temperature of the fuel pins in the inner and outer rings 
of the high power group are 1236 °C and 1177 °C respectively, which are much lower than the melting point of 
the fuel sheath. It indicates that the Canadian-SCWR can achieve “no-core-melt” concept under LOCA/LOECC. 


Keywords: Canadian-SCWR, LOCA/LOECC, No-Core-Melt, SCTRAN, Radiation heat transfer 


DOI: 10.13538/j.1001-8042/nst.26.020601 


I. INTRODUCTION 


The Canadian-SCWR is a pressure-tube type supercritical 
water reactor (SCWR), which is promising to satisfy all the 
major GIF (Generation IX International Forum) goals on en- 
hanced safety, sustainability, economics, and proliferation re- 
sistance [1]. It applies supercritical light water as coolant and 
subcritical heavy water as moderator. The high efficiency 
channel (HEC) design is employed to increase the inherent 
safety of the reactor. This helps the reactor to achieve “no- 
core-melt’” under postulated accident scenarios with a loss 
of emergency core cooling system, such as LOCA/LOECC 
(loss of coolant accident, with a coincident loss-of-emergency 
core cooling), because the radiation heat transfer inside the 
HEC and the passive heat rejection through the insulator 
into the low-temperature moderator can remove the decay 
heat. Safety analysis is required to demonstrate the feasibility 
for the continuously updated Canadian-SCWR concept. As 
LOCA/LOECC leads to the total loss of core coolant and the 
most serious accident results, its safety analysis is an impor- 
tant reference for evaluating the inherent safety of Canadian- 
SCWR. 

Simulations for evaluating and optimizing thermal perfor- 
mance of Canadian-SCWR following LOCA/LOECC was 
performed by AECL (Atomic Energy of Canada Limited) u- 
tilizing CANFLEX bundle [2]. Transient simulation was car- 
ried out with an assumed decay power variation. Effects of 
insulator properties and moderator temperature on the fuel 
cladding temperature were analyzed. The results show that 


* Supported by the Doctoral Program of Higher Education of China 
(No. 20120201110043) and Atomic Energy of Canada Limited(AECL) 
+ Corresponding author, junligou@ mail.xjtu.edu.cn 


Canadian-SCWR has the potential to significantly reduce the 
possibility of core damage frequency. Shan et al. [3] per- 
formed sub-channel analyses with ATHAS code and radi- 
ation heat transfer analyses with CATHENA code of 54- 
element Canadian-SCWR bundle. The sub-channel analysis 
results show that the maximum fuel cladding temperatures are 
761°C at BOC (Begin of Cycle) and 808 °C at EOC (End of 
Cycle). The radiation heat transfer was calculated at different 
decay power levels and the results indicate that the pressure 
tube with 54-element fuel bundle can remove about 2% of 
the rated power to the moderator through radiation heat trans- 
fer. Licht and Xu [4] provided simulation results and analy- 
ses of 78-element Canadian-SCWR bundle in the process of 
LOCA/LOECC. The results show that the temperature of fuel 
sheath with a non-porous insulator remains below the melt- 
ing temperature for less than 3% of the rated power. There 
are suggestions proposed by other authors, however, little re- 
searches about the “no-core-melt” assessment of Canadian- 
SCWR have been published. 

The previous analyses of Canadian-SCWR mainly focused 
on the radiation heat transfer, but not the effect of natu- 
ral convection during LOCA/LOECC. Meanwhile, the tran- 
sition process from supercritical pressure to subcritical pres- 
sure was not simulated in the analyses. In this paper, a two- 
dimensional heat conduction model and a radiation heat trans- 
fer model are incorporated into SCTRAN [5]. The SCTRAN 
thermal-hydraulic idealization for Canadian-SCWR based on 
the latest conceptual design utilizing a 64-element fuel bundle 
is carried out. The LOCA/LOECC is analyzed by considering 
the effects of radiation and convection in the HEC. 


Il. MODIFICATION OF SCTRAN 


SCTRAN, a safety analysis code developed at Xi’an Jiao- 
tong University, can simulate most accidents for SCWR in- 
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cluding LOCA [5], thus can be applied to simulate the acci- 
dents at both subcritical and supercritical pressures. A ho- 
mogenous model and a four-equation dynamic slip model are 
implemented into the code as optional modes for subcriti- 
cal pressure condition. Its ability to do safety analysis for 
SCWRs has been verified [5]. It has been applied to an- 
alyze the accident consequences of Chinese pressure vessel 
type concepts, such as CSR1000 [6] and CGNPC SCWR [5]. 
In order to simulate radiation heat transfer inside the HEC 
of Canadian-SCWR, a radiation heat transfer model is in- 
corporated into SCTRAN. The fuel pin absorbs different 
amounts of radiation heat from different directions, which 
results in prominent circumferential heat conduction in the 
fuel pin. Therefore, the original one-dimensional heat con- 
duction model in code SCTRAN should be updated to a two- 
dimensional one. 


A. Development of two-dimensional heat conduction model 


The differential equation of heat conduction with an inter- 
nal heat source is given as 
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where, p is density, Cp is specific heat capacity, T is temper- 
ature, t is time, À is conductivity and S is heat source, The 
first term at the right side of Eq. (1) is the energy variation 
with space while the other stands for the heat source in the 
heat structure. 

In the two dimensional polar coordinate, the conduction 
equation can be rewritten as 


OT 10 OT 1 ð OT 
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where, r is the distance in the radial direction, y is the an- 
gle in the circumferential direction . In order to reduce the 
calculation iteration and guarantee the calculation accuracy 
simultaneously, a fuel rod is divided into four sectors in the 
circumferential direction. The mesh layout for the 2D heat 
conduction is shown in Fig. 1. Points P, N, S, W and E repre- 
sent the nodes needing temperature calculation. The material 
properties and temperatures in a half mesh interval are as- 
sumed to be constants over a time interval. For the control 
volume of point P (shadowed in Fig. 1), the internal energy 
increase equals to the inner heat source plus the heat conduct- 
ed from the surrounding control volumes. Thus the discrete 
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| 
(ra + Ts) ArA@ O\ TnAn Tn 7 Tp) Tsrs5 (Tp =~) Ts) 
(C0) a Re Fe Te) = | ry, o | 3) 
_ |àe(Te — Tp) — Aw(Tp — Tw) 
| Oar. or, Ar + SAV 


where, 0 is angle in the circumferential direction. The left 
part of Eq. (3) represents the internal energy increase of the 
control volume at point P. The first term of the right side is 
the heat conducted from points N and S. The second term is 
the heat conducted from points E and W, while the last term 
denotes the inner heat source. Eq. (3) can be simplified as 


apTp = agTg + awTw + anIn + asTs + b, (4) 


where, T is the temperature for the mesh node, a and b is co- 
efficients. Every node in the mesh layout owns the same heat 
conduction equation as Eq. (4). There are five unknown node 
temperatures in each equation. Simultaneously solving the 
temperature equations of all nodes is time-consuming. For 
saving computer memory and calculation time, the conduc- 
tion equations shall be solved sector by sector. The node tem- 
peratures in one radial sector are solved together assuming 
that the node properties of the west and east sectors apply the 
values at the last iteration 


apT Ss” = ayT”) +as TE” +agT? +awTẸ T" +b. (5) 


Eq. (5) has only three unknown parameters and can be solved 
by Tridiagonal Matrix Algorithm (TDMA) with the corre- 
sponding boundary conditions [7]. When the solution of one 


radial sector is finished, physical properties of the nodes are 
updated. The radial heat conduction calculation is conducted 
in the next sector clock wise until all the node temperatures 
satisfy the convergence criterion (for each node, temperature 
difference < 1 °C between two iterations). In this way, the 2D 
heat conduction solution is transferred as many 1D solutions 
and the calculation efficiency can be greatly promoted. 


B. Development of radiation heat transfer model 


As a boundary condition for heat conduction solution, the 
radiation heat transfer inside the HEC is highly related to the 
fuel configuration. The cross section of an HEC channel is 
illustrated in Fig. 2. The radiation enclosure consists of outer 
surface of the central channel, surfaces of the fuel rods in the 
inner and outer rings, and the tube liner. The method of solv- 
ing radiation heat transfer adopted by RELAPS is referred [8]. 
The following assumptions are made: 

1. Radiation heat transfer between different elevations is 
ignored. 

2. All surfaces in the system are diffusive and grey; 

3. Radiation exchange between water steam and fuel sur- 
face do not produce a big effect on the LOCA/LOECC result, 
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Fig. 1. Mesh layout of the heat structure. 


and coolant in the pressure tube neither emits nor absorbs ra- 
diant thermal energy; 

4. Reflectance from a surface is independent of the 
reflected direction and the radiation frequency. 
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Fig. 2. (Color online) Cross-section of the high efficiency chan- 
nel(HEC). 
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1. Net radiation heat flux 


The radiosity of a surface is the radiant energy flux leav- 
ing a surface (i.e., the emitted energy flux plus the reflected 
energy flux). Energy balance for the i® surface is 


Ri =e0T +m Y Rj F Fi (6) 
=] 4 


where, R is radiation heat flux, H is Boltzmann coefficient, n 
is reflectivity, A; is area of surface i, A, is area of surface j, 
and Fj; is the view factor from surface j to surface 7. Consid- 
ering the energy conservation law, we have 


Aj Fiz = Aj Fi, (7) 


where, F;; is the view factor from surface 7 to surface j. Sub- 
stituting Eq. (7) into Eq (6), we have 


Ri = cio Tt + ni 5 Rj Fij (8) 


j=1 


The net heat flux at surface 7, Q;, is the difference between 
the radiosity of surface i and radiosity throwing to surface i 
from all surfaces, and it’s given by 


Qi= Ri- DUR Fy (9) 
j=l 


Combining Eqs. (8) and (9), the net heat flux of surface i 
is expressed by 
Qi = (Tt - Ri) (10) 
In Eq. (10), the reflectivity and emissivity are the basic 
physical properties of surface 7. By conservation of energy, 
the value sum of the reflectivity and emissivity for a sector 
equals to 1. In the calculation, the surface temperature at the 
previous time step is applied, and therefore, the radiosity is 
calculated explicitly. 


2. Radiosity solution 


On the basis of Eq. (8), the radiosity of all the surfaces can 
be calculated by: 
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The surface emissivity is regarded as constant. From [4], 
surface emissivity of the central channel is 0.34 and the emis- 
sivity of the fuel rod and the liner is 0.8. For a grey body, 
absorptivity equals emissivity and hence “n = 1 — e”, which 
would simplify the formula. Due to the central symmetry of 
the fuel assembly (Fig. 2), the pressure tube and the center 
channel are circumferentially divided into 32 sectors and each 
fuel rod is divided into 4 sectors. With the view factor of each 
surface calculated by GEOFAC code [9], the radiosity matrix 
can be solved, and the net radiation heat flux can be obtained 
using Eq. (11). 


C. Heat conduction solution considering the effect of 
radiation and convection 


The net heat flux created by radiation and convective heat 
transfer is used as the boundary condition of the 2-D heat 
conduction, 


OT 
-Apr = A(T: — Tsk) + Qi (12) 


The radiation heat flux Q; in Eq. (12) is also a function 
of surface temperature. The bulk-fluid temperature for all the 
sectors in an enclosure is assumed to be the same. The natural 
convection heat transfer coefficient inside the HEC is calcu- 
lated by Churchill-Chu correlation [10]. The forced convec- 
tion heat transfer in the HEC is mainly calculated by Dittus- 
Boelter correlation [11] at subcritical pressure and Jackson 
correlation [12] at supercritical pressure in the process of 
LOCA/LOECC. To simplify the solution, surface tempera- 
tures at the previous time step are adopted to compute the ra- 
diation heat flux explicitly and therefore the iteration solution 
of the radiation heat flux is avoided. 


I. VERIFICATION OF SCTRAN RADIATION HEAT 
TRANSFER MODEL 


Due to the lack of experimental data, verification of the 
radiation heat transfer model was done by code to code com- 
parison with Canadian system code CATHENA. As a one- 
dimensional and two-fluid thermal-hydraulic code [9], it in- 
cludes 1D and 2D heat conduction models (GENHTP). Its a- 
bility to evaluate the radiation heat transfer has been validated 
by Lei and Goodman [13]. 

The 64-element Canadian-SCWR bundle was analyzed for 
the verification. The HEC cross-section is shown in Fig. 2. 
The buddle power was 9.34 MW. The fuel pins of the in- 
ner and outer rings occupied 44% and 56% of the power, re- 
spectively. Temperature of the moderator outside the pressure 
tube was 80 °C and heat transfer coefficient between modera- 
tor and pressure tube was 1000 W/(m? K). Three steady-state 
cases with different decay heat levels (2%, 3% and 4% of the 
rated power) were simulated by SCTRAN and CATHENA, 
respectively. In the simulations, the decay heat was trans- 
ferred from the fuel rod to inner surface of the pressure tube 
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only by radiation heat transfer. Each sector in the HEC was 
marked by a number (Fig. 3). 

Table 1 lists the surface temperatures of the 64-element 
fuel bundle at 2%, 3% and 4% power levels calculated by 
SCTRAN and CATHENA. The SCTRAN results agree well 
with those of CATHENA, differing with the maximum abso- 
lute error of 6.1 °C and the biggest relative error of 0.9%. This 
shows that the 2D heat conduction and radiation heat transfer 
models developed for SCTRAN are of guaranteed calculation 
accuracy. 


Insider of the 
fuel bundle 


Outside of the 
fuel bundle 


oT 


Fig. 3. (Color online) Numbering of the HEC surfaces. 


IV. SAFETY ANALYSIS OF LOCA/LOECC 


The ability to maintain core components below the melt- 
ing temperatures in postulated accident is referred to as the 
concept of “no-core-melt’”, which is an important safety goal 
for Canadian-SCWR. The accident of LOCA with a loss of 
ECCS was simulated with the modified SCTRAN to evalu- 
ate whether Canadian-SCWR has the potential to achieve the 
goal of “no-core-melt”’. 


A. Introduction to Canadian-SCWR 


The conceptual Canadian-SCWR design possesses a mod- 
ular design that separates the coolant from the moderator, as 
in the CANDU reactors. The reactor core consists of 336 
fuel channels, each housing a 5-m long fuel assembly. It is 
designed to generate 2540 MW of thermal power or about 
1200 MW of electric power. In the conceptual design (Fig. 4), 
light water coolant enters the inlet plenum. It flows down- 
ward the central flow tube of the channel. Near the chan- 
nel bottom, the coolant exits the central flow tube, flows up- 
ward to pass through the fuel elements (fuel assembly), and 
arrives at the outlet plenum. From the outlet plenum, the high- 
temperature and high-pressure coolant is fed to the high pres- 
sure turbine. The cylindrical vessel houses the relatively low- 
pressure and low-temperature heavy water moderator. The 
main parameters of the Canadian-SCWR are listed in Table 2. 

HEC is adopted for Canadian-SCWR conceptual design. It 
consists of a pressure tube, an outer liner tube, an insulator 
and an inner liner tube. The pressure tube is surrounded by 
heavy water moderator. A potential advantage of using the 
HEC is that in the unlikely event of a LOCA/LOECC, the 
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TABLE 1. Surface temperatures (in °C) of the HEC channel calculated by SCTRAN and CATHENA at different power levels 


Power levels S1 S2 S3 S4 S5 S6 S7 S8 S9 S10 
SCTRAN 2% 1258.3 1263.8 1263.8 1236.2 1236.2 1195.8 1195.8 1125.6 1125.6 324.9 
3% 1442.8 1449.1 1449.1 1411.9 1411.9 1371.4 1371.4 1277.3 1277.3 450.1 
4% 1591.1 1598.2 1598.2 1552.3 1552.3 1512.4 1512.4 1398.9 1398.9 577.5 
CATHENA 2% 1253 1258 1258 1233 1233 1191 1191 1127 1127 328 
3% 1437 1443 1443 1408 1408 1367 1367 1279 1279 454 
4% 1585 1592 1592 1549 1549 1508 1508 1401 1401 582 


TABLE 2. Main parameters of Canadian-SCWR 


Core pressure (MPa) 25 
Thermal/Electrical power(MW) 2540/1200 
Efficiency (%) 45-50 

Inlet/outlet temperatures (°C) 350/625 

Fuel channel number 336 

Fuel bundle type 64-element 
Neutron spectrum/reactor type Thermal/pressure-tube 
Coolant light water 
Moderator heavy water 

Main coolant flowrate(kg/s) 1254 

Active core height(m) 5.0 

Cladding material Stain steel(SS310) 


Light Water 
Coolant Inlet = 


Outlet Flow « 


Moderator Outlet <= 
and Inlet = 


Fuel Assembly 
Pressure Tube 
Calandria Vessel 


Low Pressure 
Moderator 


Fig. 4. (Color online) Schematic of Conceptual Canadian SCWR 
Core. 


heat in the fuel would be transferred by thermal radiation to 
the liner tube and then conducted to the moderator through 
the insulator. The 64-element fuel bundle with two concentric 
rings is equipped inside the pressure tube. 


B. SCTRAN model 


The SCTRAN model idealization of Canadian-SCWR was 
developed based on the conceptual design. In order to simu- 
late more accurately the effect of the radial power distribu- 


tion, the 336 channels in the reactor core are divided into 
four groups (84 channels each), namely the “average power 
channels (AP)”, “high power channels (HP)’, “medium pow- 
er channels (MP)” and “low power channels (LP)” [14], as 
shown in Fig. 5 in the radial direction, In the axial direction, 
the fuel bundles and the coolant channel are divided into 10 
parts. The axial power distribution for fuel bundle refers [14]. 
Fuel pins of the inner and outer rings of each group are sim- 
ulated by independent heat structures, which exchange heat 
separately with the coolant. 


Main 
coolant 
line 
[ 
190/290 


200 


110/120 
110 


194 


Main steam 


192 i 
line 


Fig. 5. (Color online) SCTRAN idealization of Canadian-SCWR 
concept. The control volumes are numbered as 1—80. 


The heat exchange between central channel and coolant 
channel is also taken into consideration. A time dependent 
junction at 350°C and 1254kg/s, and a time dependent vol- 
ume at 625 °C and 25 MPa, are set as the boundary conditions 
of the main coolant line and main steam line, respectively. 
The moderator cooling system contains an active and a pas- 
sive moderator cooling system. In the accident conditions, 
only passive moderator cooling system (PMCS) is used [15]. 
The PMCS is simulated by volume 800 (Fig. 5), serving as 
an ultimate heat sink. It will be active automatically during 
transients. Thus, moderator temperature of 80 °C is assumed 
and a constant heat transfer coefficient of 1000 W/(m? K) be- 
tween the moderator and the pressure tube is applied [3], 
which indicates that the moderator removes the core heat 
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through natural convection. There is no active system prepar- 
ing for the Canadian-SCWR system when LOCA/LOECC 
occurs. The moderator system is assumed to keep intact in the 
whole process. Turbine stop valves are installed on the main 
steam lines, which are tripped by the “power scram” signal. 
Other initial conditions and assumptions are as follows: 

— The reactivity feedback in the core is neglected in both 
steady and accident simulations. 

— The decay heat is calculated by referring to the previous 
paper [16]. 

— The emissivity of the fuel sheath and liner tube is 0.8 
while the emissivity of central channel is 0.34. 

— The natural convection heat transfer coefficient in the 
HEC is calculated by Churchill-Chu correlation [10]; 

— Referring to the ASM specialty handbook [17], melting 
temperature of SS310 is 1400-1450°C. Thus, the melting 
temperature of the fuel cladding is set as 1400 °C. 

— Due to the radiation heat transfer effect, each fuel sheath 
had four sector temperatures. The maximum sector tempera- 
ture is applied to describe the surface temperature variation. 


C. Results for LOCA/LOECC 


The core will lose most of its inventory and suffers a cool- 
ing deteriorate condition when LOCA occurs. The heat in 
the fuel could only be removed by the passive natural convec- 
tion of the high temperature steam and the thermal radiation 
exchange between the fuel sheath and the liner tube, while 
ECCS is lost simultaneously. In this process, the radiation 
heat transfer plays a dominant role. The double-ended break 
accident at the cold leg is analyzed because it is the most se- 
vere LOCA accident. As the HP group is of the largest power 
fraction and the highest cladding temperature, it is selected to 
be the main analysis object in the current analyses. 

In the simulation, the break at the cold-leg occurred at 0.5 s. 
At 1s, the control rods started to drop and the core power de- 
creased. At the same time, the turbine stop valve was closed. 
The coolant flowrate from the main coolant pump was main- 
tained in the first 10s and then decreased to zero within 5s 
(actually decided by the pump inertia). 

Fig. 6(a) shows the coolant flowrate through the break. 
When the break occurred, large amount of coolant were 
drained out through the break at a critical velocity. The core 
pressure decreased quickly to the subcritical level because 
of the loss of coolant. The coolant flowrate in the different 
channels are shown in Fig. 6(b). The break located in the 
cold leg resulted in the reversed flow in each channel. Large 
amount of coolant of 625 °C from the outlet plenum was di- 
rected to the fuel channel, which led to a greatly degraded 
cooling condition. Therefore the fuel sheath temperatures in- 
creased quickly during this period. 

Fig. 7(a) shows the heat removed from the fuel rings and 
the heat variation in the HP group in the first 50s. Fig. 7(b) 
shows the maximum cladding surface temperature (MCST) of 
the components. At first, the nuclear power produced in the 
fuel rods cannot be totally removed to coolant and the MCST 
of both rings increased sharply. With the power decreasing 


Nucl. Sci. Tech. 26, 020601 (2015) 


30 


8 
(a) 
7@ 
D 
T 2 "S 
£ 58 
o => = 
5 4% 
8 2 
Š 10 3§ 
g ma Š 
[=] 2 č 
© € 
15 
[e] 
‘ 06 
0 10 20 30 40 50 


Time (s) 


— LP channel 


—— AP channel 
— MP channel 
—— HP channel 


Coolant flowrate (x 100 kg/s) 


0 10 20 30 40 50 
Time (s) 


Fig. 6. (Color online) The core pressure and break coolant flowrate 
(a), and coolant flowrate of the four flow channels (b), in the first 
50s of LOCA. 


and the reversed flow being established, MCSTs kept a slow- 
er increasing trend, but dropped a little at about 15s. How- 
ever, when the coolant in outlet plenum was exhausted, the 
cooling conditions in the coolant channel deteriorated. The 
channels were filled with high-temperature steam and the ra- 
diation heat transfer became dominant in the HEC at around 
30s. The heat transferred by radiation and natural convection 
at this time was not large enough to remove the decay heat of 
the inner rings, while the decay heat of the outer rings can be 
totally removed. MCST of the inner ring started to rise again 
at around 15s while that of the outer ring kept decreasing. 
Due to protection from insulator, the pressure tube surface 
remained at about 200 °C. 

With time went on, the fuel cladding temperature of the 
inner ring increased continuously while the outer ring expe- 
rienced an MCST decrease. As shown in the Fig. 8, MCST 
of the inner ring of the HP group started to drop slowly af- 
ter reaching the peak value of 1236°C. The fuel rod would 
reach a steady state when the heat removed by radiation and 
convection equals to the decay heat. However, the decay heat 
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Fig. 7. (Color online) Power decay and heat removed from the fuel rings (a) and MCST in the HP group (b), in the first 50 s of LOCA. 


decreased slowly with time, which resulted in the continuous 
decrease of the MCST. Likewise, the behaviors of the other 
three power groups are similar to that of the HP group, which 
is of the highest MCST because of its highest power fraction. 
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Fig. 8. (Color online) MCST vs time in all the coolant channels. 


According to the above analysis, the decay heat of the 
Canadian-SCWR fuel pins in the process of LOCA/LOECC 
can be transferred to the moderator outside the pressure tube 
thourgh natural convection and radiation heat transfer. The 
fuel rods of the two rings experience different sheath temper- 
ature variation. The fuel rods in the outer ring is closer to the 
pressure tube and its radiation heat exchange with the pres- 
sure tube is enhanced. While the radiation heat transfer for 
the fuel rods in the inner ring is worse than that of the outer 
ring. That is the reason why fuel rods in the inner ring ex- 
perience a longer period of sheath temperature increase. The 
peak cladding temperature of the two rings are 1236.3 °C and 
1177.3 °C respectively, which are much lower than the melt- 
ing temperature of the fuel cladding. 


[1] Leung L K H. Termalhydraulics and safety programs in support 
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V. CONCLUSION 


This paper focuses on the assessment of the inherent 
safety of 64-element Canadian-SCWR. A two-dimensional 
heat conduction model and a radiation heat transfer model 
are developed and incorporated successfully into SCTRAN. 
Meanwhile, the validity of the newly developed models for 
SCTRAN is verified by CATAENA, which is a system code 
developed by AECL. 


A SCTRAN idealization of the Canadian-SCWR concep- 
tual design has been developed. Using the SCTRAN ideal- 
ization developed for Canadian-SCWR, a simulation of the 
steady-state thermal-hydraulic conditions is performed. In 
LOCA/LOECC, with total loss of core coolant inventory, ra- 
diation heat exchange among sheaths and the natural con- 
vection of the high-temperature steam can removed certain 
amount of core decay heat. On the basis of conservative as- 
sumption, the maximum fuel sheath temperatures in the inner 
and outer rings of the HP group are 1236°C and 1177 °C, 
respectively, which are lower than the melting point of the 
cladding material. The results show that the Canadian SCWR 
is capable of achieving the design object of “no-core-melt” 
under LOCA/LOECC. 
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